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We consider a model of two (fully) compact polymer chains, coupled 

through an attractive interaction. These compact chains are represented by 
Hamiltonian paths (HP), and the coupling favors the existence of common 
bonds between the chains. Using a (n = component) spin representation 
for these paths, we show the existence of a phase transition for strong cou- 
pling (i.e. at low temperature) towards a "frozen" phase where one chain is 
completely adsorbed onto the other. By performing a Legendre transform, 
we obtain the probability distribution of overlaps. The fraction of common 
bonds between two HP, i.e. their overlap q, has both lower (q m ) and upper 
(qm) bounds. This means in particuliar that two HP with overlap greater 
than qu coincide. These results may be of interest in (bio)polymers and in 
optimization problems. 
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I. INTRODUCTION 



The study of two coupled identical systems is a familiar topic in Statistical Physics, both 



exploration of phase space. In this note, we wish to extend these studies to the case of two 
compact polymer chains. The coupling between the chains is chosen as an attractive term 
which favors the existence of common monomers (bonds); in other words, we consider the 
"adsorption" of a chain onto the other. By a Legendre transform, this model maps onto the 
overlap probability distribution of the chains. 

For simplicity, we will consider fully compact chains. This allows us to use a simple 
spin representation || to describe the chain properties through Hamiltonian paths. Beside 
its interest in homopolymer physics, this model may yield some insight in the sequence- 
stucture relationship in proteins, or in some related optimization problems such as the 
travelling salesman. The outline of the paper is as follows: we briefly recall in section [IT] 
the connection between a single polymer chain and an n-component spin system, in the 
limit n = 0. The extension to the problem of two interacting chains (section |TXTf) is then 
straightforward, and will be solved at a mean field level in section [IV]. We finally mention 
some possible consequences of the model. 



The connection between polymer physics and spin systems may be presented as follows 
f|. Let us consider, on each site r of a d- dimensional cubic lattice, a spin variable S r , with 
n components. The normalization is such that S% = n. Consider the sum: 



where the variable K denotes the fugacity of a monomer, and the product runs over neigh- 
bouring pairs of sites (i.e. bonds) of the lattice. In eq ([!]), dfi is the normalized integration 
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in equilibrium [|I] and non equilibrium contexts [0] . It is commonly considered in the physics 
of glassy systems where it is a substitute to the time honored conjugate field in the 



II. THE SINGLE CHAIN PROBLEM 




(1) 



measure on the (n — 1) dimensional sphere of radius \fn. 
Due to the normalization of S r , it is easy to see that 

J dfx(S r ) ex V (H r S r ) = 1 + S- + 0(n) (2) 
so that, by taking derivative with respect to (a = 1, ■■■,n), we obtain: 

J dn(S r ) 1 = 1; / dfi(S r ) (S?) 2 = 1 
All higher powers of S r are at least of order n. 

Expanding eq ([!]) in powers of K, and using the previous remark, we see that Z(K) can 
be viewed as a sum over all closed loops, each closed loop of length I contributing a weight 
nK l . In the limit n — > 0, only the single connected loops survive, hence the possibility to 
represent self-avoiding walks (SAW) on the lattice by eq ([!]). A very useful rewriting of eq 
(B) ^ 

Z(K) = J J] d^(S r ) e^ E <^> KSrA ** § r' (3) 

r 

where the operator is a lattice A function (A^ = 1, if r and t > are neighbouring sites, 
and otherwise). Note that the sum in the exponential term of eq ([3]) is over the sites of the 
lattice. Applying the familiar Stratonovich-Hubbard transformation to eq @, we introduce 
n = component fields r and get 

Z{K) = Ji[d$ r e" 1 E (-') $rA ^ r ' J J] dfi(S r ) e^^l (4) 

f r 

Using eq @, we get 

z(k) = J n#; e-^)* rA ^> jj h + (5) 

r f 

So far we have not specified the spatial extension of the SAW: Z(K) is the grand partition 
function of the chain, so that the number M of sites of the lattice is not related to the 
number N of bonds of the SAW. We now require that fully compact configurations are 
the only configurations present in eq (||), i.e. M = N. This full compactness requirement 
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amounts to keep the term proportional to nK N in eq ([5]). The coefficient of this term is 
simply the number M of self avoiding fully compact configurations (also called Hamiltonian 
paths or HP), that is 

n = k-« (Mi)^ = i im „.„i/ri#; e-*^^'^ n (f ) <« 

Note that M is simply the canonical partition function for a HP of N bonds. Performing a 
homogeneous saddle point approximation on in eq (§) yields 

tf=(f)" (7) 

where 2 = 2d is the coordination number of the lattice. Note that in the context of Hamil- 
tonian paths, an homogeneous solution implies that one deals with periodic boundary con- 
ditions [|J. 

Equation (P) holds for an arbitrarily connected graph with adjacence matrix A^/, but 
the saddle point evaluation is a priori valid only for large enough z. Indeed, for the fully 
connected graph [z = N), the number of HP can be directly estimated (A/ - = N\/(2N)), 
and eq (0) reduces to the Stirling formula. 

III. THE COUPLED CHAINS PROBLEM 

We now extend the above approach to the case of two interacting self avoiding chains. 
By "interacting chains " , we mean that the presence of common (or doubly occupied) bonds 
between the chains is thermodynamically favored, and we are interested in counting the 
number of configurations of two such HP. 

Following eq ([!]), we denote by K the fugacity of a monomer and consider the sum 

Z 2 (K, A) = f Y[dfi(S r )dfi(a r ) ]J (l + K(S r S r > + oyoy) + \K 2 (S r S r/ )(a r a r ,)) (8) 

where a r has the same properties as S r (a 2 = n), and the parameter A is a measure of the 
interaction between the chains (see below). The contribution of doubly occupied bonds to 
eq (|D is proportional to A^ 12 , where Nu is the number of these bonds (implying A > 1). 



Following exactly the same steps as in section ||, we have 

Z 2 (K, A) = J IJ d^S r )dfi(a r ) e^ E (^) A ^°^ (9) 

r 

with 

G 2 (S, a) = K(S r S r > + a r a r ,) + (A - l)K 2 (S r S r >)(a r a r ,) (10) 
Performing the Stratonovitch-Hubbard transformation on eqs.@ and ([TOD yields 

Z 2 (#, A) = / 1] #r#, II <kaf>(r) e"^ A ^ | J] dfi(S r )dfi(a r ) e^Y*** (11) 



with 



a/3 



Aj*j* = 4>r<Pr> + IprA' + /2 Qa^Qapi^) ( 12 ) 

a/3 



and 



B r = S r $ r + ff r f r + ^v 7 !^! ]T c? a/3 (r ) S?a? (13) 

a/3 

— # — * 

Note that the previous transformations require two n component fields <p r and if) r , and an 
n x n matrix q a p{r) , with n = 0. One may now perform the integration over the original 
spin variables S r and a r in eq (|TT|) . We obtain 

Z 2 (#, A) = J ]Jd$ r dAl[dq a p(r) e ^ E (^) A ^ JJCr (14) 



a/3 



where 



Cr = i + + + + ^v^EC^W^f + ^(a - 1) E^W (is) 

a/3 a/3 

So far, we have not specified the spatial extension of the chains. The sum Z 2 (K, A) is grand 
canonical with respect to both singly and doubly occupied bonds. Following the previous 
section, we now enforce the full compactness constraint for both chains, by keeping the term 
proportional to n 2 K 2N in eq (13). Defining by Z(X), the total number of HP of the two 



interacting chains, we have 



Z{\) = K~ 2N (^^] (16) 

Technically, the extraction of the term proportional to K 2N in Z 2 {K, A), amounts to keep 
on each site r the terms proportional to K 2 in eqs.(|l4]) and (p~5f ) . We may therefore write 



Z(X) = lim^o-^ J J] dlA J] dg Q/3 (r) e _l TJ P r (17) 



jL f T\ AX Ant J]V~ J ^ ° ^ Efr,r'-| ^^rP'- 4 ??' 
a/3 

where is given in eq ([O) and 



v r = ^f + v^E^M^f + ^(A - i) £<&(r) (is) 

a/3 ~" a/3 

Note that 2(A) is still grand canonical with respect to the number N 12 of common bonds 
between the two HP, since, as mentionned above, the term with iV 12 common bonds in Z(X) 
yields a contribution proportional to A^ 12 . 

One may also interpret eq flTB]) as giving the partition function of two HP of N bonds, 
with an attractive interaction energy e favoring common bonds, at temperature T, that is 

Z(e)= e ^ Nl2 ( 19 ) 

(HPi,HP 2 ) 

with (3 = if. The identity between eqs.(|l6|) and (|19|) , i.e between Z(\) and Z(e), yields the 
familiar result A = e@ e . Moreover, equation (|19"D allows us to derive bounds for the partition 
function Z(e) (or Z(X)), namely 

Me PeN < Z{e) < Af 2 (20) 
where TV is given in eq (|7|). From now on, we will set e — 1. 

IV. SADDLE POINT APPROXIMATION AND PHASE TRANSITIONS 

Since an exact evaluation of eqs (|T7| , |T8| ) seems to be out of reach, we will use a saddle 
point approximation with respect to the variables <f) r , ip r , and q a p(r). This saddle point will 
be further restricted to be space independent (0 r = 0, ip r = if), q a p{r) = q a /3), again implying 
periodic boundary conditions. Setting Z(X) = e~ Nu) W ? we get 



u(X) = Min ( ^ >qa0) ( -L$» + ^ + qlp) - LogX? ] (21) 



a/3 



where 



0V 



V = Zf- + v^Efg^/ + - (A - 1) £ (22) 

a/3 a/3 

The saddle point equations are easily solved, by introducing the quantities $ = <p 2 , \I/ = 
■?/> 2 , R = (j> a Qa/3^ \ and Q = 2 a( g Denoting the saddle point values with a subscript 
0, we obtain 

$ _ |$2 + VA^TRo 



2d £> n 



Ro _ (A - 1)R + V\~=T& 

2d Pn 



Qo = (A - l)Qp + y/A^TRo 
2d V 

together with $o = ^o- in eqs. ( p3| - |25|) , the denominator is given by 



(23) 



(24) 



(25) 



One then gets 



Do = M + ^TR„ + ^-y^Qo (26) 



At this point, one has to resort to a numerical solution of the saddle point equations. 
Note that the possible solutions obey simple equalities or inequalities, such as: $o, Qo > 
0, $ + Qo = 4d, R 2 = Q $g. 

Our results can be interpreted in two ways. They first describe the (thermal) properties 
of two HP coupled via eq flUf): LogZ(A) = —Nuj(X) is then, up to a temperature factor, the 
free energy (with A = e 13 ). On the other hand, a Legendre transformation with respect to 
LogA, gives information on the number Nu of common bonds between the HP, i.e. on their 
overlap properties. For convenience, we define = Nq, and 



(28) 



A saddle point evaluation of eq fl2~8]) gives 



dcu(\) 
dLogX 



(29) 



Using eqs (|23| - p7D , we obtain 



A Qo 
A - 1 4d 



(30) 



We now present our results along both ways. 

A. Thermal properties of the coupled Hamiltonian paths 

Unless otherwise specified, our results are given for d = 3, and the main parameter of 
this section is the temperature T (T = j^j). The saddle point equations ( |2~3| - |27|) yield the 
free energy fo(T) = Tuq^t) as a function of T (see Figure [I]). 




FIG. 1. The free energy / as a function of temperature. The saddle point result /o (full line) 
crosses the free energy 7 of the fully adsorbed phase (hatched line) at point A (T = T3). 



Solving numerically the above equations, we get the following results 
1) T x = 00 



S 



In this case, we find a finite fraction of common bonds q$ — h — ~. This value 
corresponds to a random choice of a common bond among the d = 3 lattice dimensions, and 
is the smallest possible value q m of the overlap. 

2 ) T * = l^+i) - °' 7213 

One then has $o = Qo = 2d = 6. This point is in some sense a disorder point || where 
the values of $o (linked to the entropy of a single HP) and Qo (linked to the overlap of the 
two HP) cross. It corresponds to an overlap = 4^- = |. 

3) T 3 ~ 0.5846: the complete adsorption transition. 

At this point, the two HP system "freezes" into a single HP: one may also say that one 
chain gets fully adsorbed onto the other. Of course, the entropy does not vanish below T 3 , 
but becomes equal to the entropy of a single HP, see eq (0). Note that this solution is not a 
saddle point solution. It must nevertheless be taken into account, since its free energy per 
monomer which reads 

T (T) = -l-TLog(~) (31) 

is an upper bound (see eq (^)). The transition temperature T 3 is thus defined by 7(T 3 ) = 
foiTz) (see Figure |l|). The mecanism of this first order freezing transition is analogous to 



the one studied in ref. flO|l , in the context of polymer crystallization. Note that q jumps 
(3) 

from g ~ 0.7572 to q — 1 across the transition. From a thermodynamic point of view, we 
have thus found that the overlap fraction of two HP has a maximum value qu — 0.7572, 
beyond which the two HP coincide. 

Note that, as the space dimension d increases, the adsorption temperature T 3 decreases. 
In particular, the transition does not occur for the fully connected graph (d — > 00). 

4) In the present mean field like description, one may continue the high temperature 
branch of the free energy fo(T) beyond the full adsorption point, implying the existence of a 
metastable state (see the branch AB of Figure 0). In this respect, a particular temperature 
may be defined, where the entropy s (T) = — equals that of a single HP: 

So (T 4 ) = Log(^) (32) 
9 



leading to T4 ~ 0.4350, and in turn to q$ ~ 0.9627. Furthermore, one may argue that this 
metastable branch is defined up to the point where the (metastable) overlap parameter q 
reaches the value one. If this is correct, we find from the saddle point equation that this 
occurs for T = T 5 ~ 0.3898, with a positive entropy so(T 5 ) ~ 0.498. We do not have a clear 
understanding of the metastable branch, and more work is needed on this point. 



B. Overlaps and Legendre Transform 

As previously mentionned, our results can also be interpreted in terms of the overlap 
properties of the two HP. A convenient function to characterize these properties is the 
overlap probability distribution V(q) defined as 

V (l) = Tfl E S(N a -Nq) (33) 

where Af is given in eq (0) and iV 12 is the number of common bonds of the HP. Eqs.(^) and 
(|33|) imply that 

V(q) = jpe(q) (34) 

From a strictly thermodynamic point of view, we have the result that the overlap prob- 
ability distribution V(q) is defined only for q m < q < q M , and for q — 1. The existence of 
a metastable branch beyond qM is not easy to interpret: the saddle point evaluation (|29|) 
of the Legendre transform (|28"D then becomes ill-defined, since \(q) becomes a multivalued 
function. 

For q m < q < qM, the relation between A (i.e. temperature) and q (i.e. overlap) can be 

inverted through eq (^). There is thus a one to one correspondance between overlap and 

1 

temperature in the region 1 < A < A3 = e T 3 ~ 5.532. The thermal properties of the coupled 
system imply that, for 1 < A < A 3 , we have V(q) = 6(q — q), where q is given by eq (|30|). 
On the other hand, for A > A3, we have V(q) = 5(1 — q). 

Another quantity of interest is the entropy, considered as a function of the overlap q. It 
is given by 
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s(q) = 1 Log6(g) (35) 
Performing the saddle point evaluation on A in eq fl28|), we get 

s (q) = -o; (A) - gLogA (36) 



where A = X(q) is given by eq (|30|). The phase transition for q = qu can be interpreted 
(see Figure 0) as a Maxwell construction since the results of the previous section can be 
rewritten as 

So(<?m) - «o(l) = (1 - 9m) LogA 3 (37) 



1.6 
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FIG. 2. The entropy s(q) as a function of the overlap q. The full line denotes the saddle point 

result so(q); the hatched line is the result of the Maxwell construction with the fully adsorbed 
state. Note that AB corresponds to the metastable branch of Figure |l|. 



Beyond the homogeneous saddle point approximation, our results in the region qn < q < 
1 are coherent with a phase coexistence picture between a fraction x of the phase q = qu 
and a fraction 1 — x of the fully adsorbed phase q — 1. Such a phase coexistence will give 
rise to an effective overlap qu < q e ff — x Qm + (1 — x) < 1. 
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V. CONCLUSION 



We have considered the overlap and/or thermal properties of two coupled Hamiltonian 
paths (HP), in a homogenous saddle point approximation. We have found (for d = 3) a 
phase transition at finite coupling, between an entropy dominated phase and a completely 
adsorbed phase, where the two HP have overlap q — 1. This phase transition disappears 
when d —>■ oo. Our results may be relevant in a proteic context (NMR, structure alignement). 
If the number of constraints exceeds a certain threshold, the existence of qu suggests that 
a single spatial structure may survive. For longer polymers, the existence of a first order 
transition, raises many questions (metastability, glass transition,..). Finally, we remark that 
the coupled chains system of this paper undergoes a phase transition at a finite value of the 
coupling constant A, whereas in spin glasses, the two replica system undergoes a (spin glass) 
transition when the coupling constant vanishes @^]. 

S. Franz thanks the Service de Physique Theorique, Saclay, for its kind hospitality during 
the elaboration of this work. 
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